Nongaussianity from Tachyonic Preheating in Hybrid Inflation 
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In a previous work we showed that large nongaussianities and nonscale-invariant distortions in 
the CMB power spectrum can be generated in hybrid inflation models, due to the contributions of 
the tachyon (waterfall) field to the second order curvature perturbation. Here we clarify, correct, 
and extend those results. We show that large nongaussianity occurs only when the tachyon remains 
light throughout inflation, whereas n — 4 contamination to the spectrum is the dominant effect 
when the tachyon is heavy. We find constraints on the parameters of warped-throat brane-antibrane 
inflation from nongaussianity. For F-term and D-term inflation models from supergravity, we obtain 
nontrivial constraints from the spectral distortion effect. We also establish that our analysis applies 
to complex tachyon fields. 

PACS numbers: 11.25.Wx, 98.80.Cq 



I. INTRODUCTION 

In the simplest models of inflation, the primordial den- 
sity perturbations have a negligible degree of nongaus- 
sianity. The parameter /jvl which characterizes non- 
gaussianity is of the order of \n — 1| <C 1 (where n is 
the spectral index) in conventional inflation models 
4], whereas the current experimental limit is |/jvl| r$ 100 
5j; one can additionally characterize the nongaussianity 
using the trispectrum, which is also small in conventional 
models @. Nevertheless, there have been intense theo- 
retical efforts to find models which predict observably 
large levels Q (see Q for a review). It has been difficult 
to find examples which give large /nl- In single field in- 
flation models a small inflaton sound speed is necessary 
to achieve large nongaussianity [9(, as in the models of 
1C , [llj , unless the inflaton potential has a sharp feature 

a. 

The simplest multi-field models do not seem to give 
large nongaussianity [HI , though it is not clear if this is 
true also of more complicated models. Thus it is quite 
significant that one of the most prevalent classes of mod- 
els, hybrid inflation [l4j, is able to yield large nongaus- 
sianity for certain ranges of parameters [1 51 ] . The ef- 
fect is due to the growth of the waterfall field — tachyonic 
preheating — which contributes to the curvature pertur- 
bation, and hence the temperature anisotropy, only start- 
ing at second order in cosmological perturbation the- 
ory. We show that, depending on the values of certain 
model parameters, two interesting effects are possible: 
n = 4 distortion of the spectrum or large nongaussianity. 
The same effect was observed in [l||, though tachyonic 
preheating after hybrid inflation was not considered in 
that paper. 1 Nongaussianity from preheating has also 



been studied in [17|-[22]. The calculations are compli- 
cated, and required numerical integrations over time and 
wavenumbers; hence the results are not immediately in- 
tuitive. One of our goals in the present paper is to 
give a better understanding of this novel effect, and to 
present some new results concerning the application of 
these results to popular models of hybrid inflation includ- 
ing brane inflation [23| and P-term inflation (24j], which 
is a synthesis of supergravity inflationary models inter- 
polating between F-term and D-term inflation. 



We begin by reviewing the results of flj| in section 
HI1 In section IIIII we apply these results by establishing 
constraints on the parameters of hybrid inflation, com- 
ing either from the production of large nongaussianity, 
or from nonscale-invariant contributions to the spectrum 
(as opposed to bispectrum). These results extend and 
correct our previous limits [151 ]. We then adapt them to 
the cases of brane-antibrane inflation in section IIVI and 
P-term inflation in section |Vj We further extend our 
analysis to the more realistic case of a complex tachyon 
field in section [VTl showing that the extra components of 
the tachyon add in a simple way and amplify the real- 
field results by factors of order unity. Conclusions are 
given in section IVTT1 Appendix A gives details about the 
matching between early- and late-time WKB solutions of 
the tachyon fluctuation mode functions, while appendix 
B gives details about the source term of the curvature 
perturbation for complex tachyons. 



1 Reference 16| studied the behaviour of metric perturbations dur- 
ing preheating in the model V = A</> 4 /4 + # 2 2 x 2 /2-|-A'x 4 /4 and 
found that an n = 4 contamination of the spectrum is generated 



when the x field is heavy throughout inflation while large non- 
gaussianity is possible when \ field is light. In this paper we 
consider a different model, tachyonic preheating, finding results 
which are qualitatively similar. 
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II. REVIEW OF PREVIOUS RESULTS 



A. Hybrid inflation 



panded metric and Einstein equations can be found in 
[15| . The matter content of the theory is expanded in 
perturbation theory as 



The hybrid inflation model which we study is defined 
by the potential 



V(<p,<r) = j 



'a 2 - v 2 ) 



-% 2 



9 2 2 

— tp a 
2 * 



(1) 



where tp is the inflaton and a is the tachyonic field. Its 
mass depends on tp as m 2 = — Xv 2 +g 2 <p 2 , which changes 
sign when <p reaches the critical value ip c = (y/X/g)v. 
At this time, fluctuations in the tachyon field start to 
grow exponentially. This phase of exponential growth is 
called tachyonic preheating [H[-[28j]; see also [2!| for a 
discussion of the general theory of preheating and [3(| 
for a different type of tachyonic preheating. 

During the early stages of preheating, before the fluc- 
tuations have become nonperturbatively large and before 
the backreaction has set in, the expansion of the universe 
will still be approximately de Sitter. Once the tachyon 
fluctuations become sufficiently large their backreaction 
modifies the expansion of the universe and brings infla- 
tion to an end. This happens at a time N* = Ht* when 
the fluctuations in a grow to a certain value, 

d 3 k 



(5a 2 (A*)) 



(2tt): 



(2) 



N=N+ 



where is the mode function for the fluctuations (dis- 
cussed below). This happens at some time after the onset 
of the instability. For a wide range of parameters (in- 
cluding the values originally considered in [l4j]) one has 
JV 4 < 1 so that the symmetry breaking completes on a 
time scale short compared to the Hubble time. This is 
the usual waterfall condition regime of hybrid inflation. 
In the present work we will consider both the possibilities 
that N± < I and also AT* > 1. 

We find it convenient to measure time in terms of num- 
ber of e-foldings, taking N = to coincide with the onset 
of the instability, when m 2 — 0, and A* to be the end 
of inflation, defined by Horizon crossing occurs at 
some Ni < 0, so the number of e-foldings since horizon 
crossing is N e = AT* — Ni. We determine N e using the 
standard relation 

' (3) 



N e = 62 



, /I0 16 GeV\ I 
In r-7-j — In 

V F 1/4 / 3 



Pr.h. 



with the energy density at reheating (p r .h. ~ T* h ) as- 
sumed to be limited by the gravitino bound T r h. £ 10 10 
GeV, though we have checked that this assumption has 
little effect on our results. Given A 7 * and N e , Ni is de- 
termined by Ni = N* — N e . 



B. Second order curvature perturbation 

We work up to second order in perturbation theory, 
employing the longitudinal gauge throughout. The ex- 



ip(r,x) = po(r) + 5^p(r,x) + -5^p(t,x) 

o-(t,x) = J (1 V(r,x) + i<5 (2 V(T,x) 

As discussed in fl5j , we are justified in dropping the ho- 
mogeneous backgound of the tachyon field (o-(t,x)) = 
o~o( T ) = 0. Conformal time, r, is related to cosmic time 
as dt = adr. We denote derivaties with respect to con- 
formal time as /' = d T f and with respect to cosmic time 
as / = d t f. 

Similarly the gauge invariant curvature perturbation, 
C, is expanded in perturbation theory as 

c = c (1) + \t {2) 

Because o~o = the first order contribution f^ 1 ) is iden- 
tical to the standard result from single field models. We 
split the second order curvature perturbation into a com- 
ponent which is due to the inflaton field and a component 
which is due to the tachyon field as 

= + C (2) 

(2) 

The second order inflaton curvature perturbation, Q, , 
coincides with the in single field models. This con- 
tribution has been previously computed and is known to 
be small and conserved on large scales 

(2) 

The quantity of interest is Q- , the tachyon curvature 
perturbation. Beyond linear order in perturbation theory 
there are nonadiabatic pressures in the model which will 
source the time evolution of on large scales. The 

(2) 

contribution Q- is the term which is amplified during 
the preheating phase and which will come to dominate 
£( 2 ) at late times. We therefore focus on \ since any 
significant nongaussianity will arise due to this term. 

One of the principal results of [lj| was the computation 
of the tachyonic contribution to the second order tachyon 
curvature perturbation in terms of the first order tachyon 
fluctutations S^'a: 



'(2) 



K 

e 



dr' 



(5^0-' j 



i\2 



H{rf 



(4) 



where k 2 — M^ 2 — 8itGn, e is the slow roll param- 
eter, e = \M 2 (V'(p)/V) 2 , t is the conformal time, 
t = — [Ha(l — e)] _1 , H is the conformal time Hubble 
parameter, H = 1[t(1 — and all factors in the in- 

tegrand are evaluted at r' unless otherwise indicated. In 
deriving ([5]), we have performed partial integrations in 
which surface terms at the initial time were dropped; 
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hence ((4]) is only valid for perturbations which are dom- 
inated by the tachyonic growth at late times. For such 
perturbations, there is little sensitivity to the value taken 
for Tj. An analogous result was derived for hybrid infla- 
tion (not considering preheating) in [3l| . Metric per- 
turbations during preheating have also been discussed in 

(2) 

Since Q depends only on the first order tachyon fluc- 
tuation 5Wcr, and not on S^a 2 we can drop the super- 
script (1) and denote S^a by 5a in subsequent text. 

The expression satisfies an important consistency 
check, namely that it is a local expression. Nonlocal op- 
erators A~™, powers of the inverse Laplacian, arise at in- 
termediate steps in the calculation, using the generalized 
longitudinal gauge, which separates metric perturbations 
into scalar, vector and tensor components. In the process 
of decoupling these to solve for the curvature perturba- 
tion, one must apply A -1 . The nonlocal terms should 
cancel out of physical quantities, similarly to electrody- 
namics in Coulomb gauge. The second order curvature 
perturbation is related to the observable CMB temper- 
ature fluctuations in a nontrivial way, so this by itself 
does not prove that £( 2 ) must be local. However [36[ has 
recently shown that under the conditions present in our 
model, C'" 2 '' should indeed by local. 

Using it is possible to compute tachyonic contri- 
butions to the spectrum and bispectrum of the curvature 
perturbation, 



(2) A2) \ _ 



>(2) A2) 
^>fci S& 2 "»fc 3 



= 6(ki + fe) S(h) 
_ (5( 3 )(fci + fc 2 + fc 3 ) 

(27T)3/2 



(5) 



B(ki, k 2 , k 3 ) 



(6) 



C. Tachyon mode functions 

To compute the correlations in (©-j6|), we express the 
tachyon fluctuation in terms of creation and annihilation 
operators, 



Sa(x, N) 



d 3 k 
(2tt) 3 / 2 



a k £k(N)e 



ik-x 



h.c. 



(7) 



where the mode functions obey the linearized tachyon 
equation of motion. To make this equation more 
tractable, we have approximated the tachyon mass de- 
pendence on time (N) as being linear, m 2 = —cH 2 N, so 
that 



d2 r o d x 
_Sa k + 3—Sa k + 



k 2 e~ 2N 



cN 



Sa k = 



(8) 



2 Indeed, as was shown in |l5j . S^cr decouples from the gauge 
invariant quantity up to second order in perturbation theory. 



where k = k/H. We found that this technical assump- 
tion is nearly always satisfied for model parameter values 
consistent with the near scale- invariance of the CMB fluc- 
tuations; furthermore it is always satisfied for parameters 
which lead to large nongaussianity. The term cN is pro- 
portional to the first term in the Taylor series for e 2ri N — 1 , 
where rj = 4M 2 m 2 i /(Aw 4 ) is the usual slow roll parame- 
ter for ip. We thus demand that 2r;|iV| <C 1 throughout 
inflation. Notice that inflation is ended by the tachy- 
onic instability, not by the failure of the slow roll condi- 
tions. The coefficient c is given by c = 2i]\v 2 / H 2 , and 
H 2 = y/(3M 2 ), with V = j\v 4 . Using the COBE nor- 
malization V/(Mp€) — 6 x 10~ 7 to eliminate the inflaton 
mass ro», we find that c = 2.2 x 10 4 g M p /v. 

Although eq. {8]| has exact solutions in terms of Airy 
functions when k — 0, for general k no closed-form so- 
lutions exist. We therefore constructed solutions, using 
the WKB approximation, or alternatively the adiabatic 
approximation. The WKB approximation is valid in the 
limit of N — > ±oo; these solutions are matched to each 
other at N k , a fc-dependent intermediate value of N. 
They have the form 



/ „ \ -1/2 / 

\2Hk 3 ) (l + ike- N 



N <N k 



b k e~i N +^ z / (1 + H)- 1/4 , N>N k 



(9) 



with 



b k = 



i-t^NT\ (1 + N) 1 / 4 



/ 2i?(c|^|) 3 / 4 exp(^; 



3/2 



(10) 



where z = (l + |ciV) , z k = (l + fcAj.) , and the dividing 
time between the small- and large- N behavior for a given 
mode is implicitly defined by N k = \n(k/y/c) — In y/\N k \- 
The matching time N k is discussed in some detail in ap- 
pendix A. The alternative method, using the adiabatic 
approximation, will be reviewed in section See also 
[37| for a discussion of the solutions of the mode function 
equation. 

The solutions of {H have been discussed in great detail 
in [l5j . However, a few comments are in order about the 
solutions ©. At early times N < N k the gradient term 
in the Klein-Gordon equation dominates over the mass 
term, k 2 /a 2 > |m 2 |, and the resulting mode functions 
look just like the solutions for a massless field in de Sit- 
ter space. These ultraviolet modes are redshifted by the 
expansion of the universe into the instability band where 
the mass term dominates the dynamics |?n 2 | > k 2 /a 2 . 
(Because of the time dependence of m a the modes may 
reenter the massless regime for a brief period of time; we 
have verified that this does not alter any of our results.) 
In this infrared regime (where the mass term dominates) 
the Airy function solutions are appropriate. We have 
checked the solutions ([9|) against numerical solutions of 
©, and found good agreement. 
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D. Integrated results 

Using the solution (J9J) in (4]), and going to the limit of 

(2) 

vanishing wave numbers, Q takes the form 
^ 2 ' d 3 p 



(apbp + a\ b- p ) 



e J (2tt) 3 / 2 
x / f(c,N,N*)dN (11) 

Jma,x(N p ,Ni) 

where f(c, N, N*) is given by 



f(c,N,N*) = e 

9 
4 



(l - e 3 ^- w *)) 

cNe HN-N,) 



r z -\- 



2c sign(z) 



27(1 + \z\ 



(12) 



It should be noted that the dominant time-dependence is 
determined by the combination — 3iV+^z 3 / 2 in the expo- 
nent. The e~ 3N decay factor is typical of massive modes, 
which redshift as 5ak ~ a -3 / 2 . The e^ z3/ growth factor 
is a result of the tachyonic instability. As we noted in the 
discussion following eq. (H|), (fTTj) is valid only when the 
late-time behavior dominates. Therefore an important 
consistency condition for all of our analysis is that 



2c Z * 



3/2 _ 



3/2 



> 3|JVi| 



(13) 



If this condition is violated then the preheating is not 
playing any role in the dynamics and we can safely as- 
sume that no significant nongaussianity is produced. 

Using (fT2"|) . the correlation functions which give the 
spectrum and bispectrum from ([5]) and © can then be 
computed as 



and 



(2tt) 3 



max(JVp,JVi) 



dN f(c,N,N») 



dNf(c,N,N*) 

max(JV p ,iVi) 



(14) 



(15) 



The tachyonic contribution to the spectrum cannot 
exceed the experimentally inferred inflaton power spec- 
trum, y/P<p(k) = 2tt x lCT 5 *;- 3 / 2 , leading to the 
bound S < P v . Moreover the bispectrum is re- 
lated to the nonlinearity parameter f^L by B = 
— ^Jnl (Pip(ki) P v (fo) + perms), which at equal mo- 
menta ki leads to 



Tnl 



5 B 
18 Pi 



(16) 



The current experimental constraint is |/atl| ~ 100. By 
analogy, we also define a parameter fa for the spectrum 
as 



f S 



(17) 



and demand that < 1. 



CONSTRAINTS ON HYBRID INFLATION 
PARAMETER SPACE 



III. 



By numerically evaluating the integrals (fT4")) and (fTJ 
and applying the experimental limits on the inflaton 
power spectrum and bispectrum, we were able to find 
constraints on the hybrid inflation parameter space. We 
update these bounds in the present section. 



A. The issue of scale invar iance 

In [ll| we noted that it is possible for the tachyonic 
contributions to the spectrum or bispectrum either to be 
nearly scale invariant (S ~ 1/fc 3 , B ~ 1/fc 6 ), or else to 
badly violate scale invariance (S, B ~ k°). In the lat- 
ter case, the spectral index for the tachyon contribution 
to the two-point function is n = 4. However we did not 
clearly differentiate between these two regimes in the lim- 
its presented in [1 51 ] . an omission which we rectify here. 

The two regimes, scale-invariant and nonscale- 
invariant, can be understood in reference to the condition 
(I13[) which must be satisfied in order for tachyonic pre- 
heating to play any significant role. There are two ways 
to satisfy (fT3j) . One is to take cTV* 3> 1, which usually 
requires c > 1. This is the regime in which the tachyon 
mass is not small compared to H during most of inflation, 
and so it corresponds to nonscale-invariant fluctuations 
of 5a. The tachyon fluctuations are Hubble damped as 
8a ~ a(i)~ 3 / 2 prior to inflation, but this suppression can 
be overcome on large scales if the amplification during 
the preheating phase is sufficiently large, which typically 
requires very small values of the self-coupling A < 1. 
This nonscale-invariant regime corresponds to a region 
of the parameter space where the waterfall condition of 
hybrid inflation is satisfied. 

The second way to satisfy (IT3l) is to take c\N\ < 1 for 
all N £ [Ni,N±]. This gives a scale- invariant spectrum 

for the tachyon and also for Q- ; , which is most easily seen 
by writing the tachyon mass-squared as \m 2 a \/H 2 — c\N\. 



3 One might consider being more conservative and imposing, say, 
< 0.01, rather than < 1, as we have done. Because 

the effect turns on exponentially fast, our exclusion plots are 
actually quite insensitive to the value assumed for fi, and /jvl- 
For example, the exclusion plots for |/jvl| < 1 are visually hard 
to distinguish from those for I/atlI < 100. 
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It is clear that if c\N\ < 1 for all N G [Ni,N*] then the 
tachyon field will have been light compared to the Hubble 
scale for all ~ 60 e-foldings of inflation which ensures a 
nearly scale-invariant spectrum for the tachyon. Also, 
in this case the instability will typically take several e- 
foldings to complete so that N* > 1. This scale-invariant 
regime corresponds to a region of the parameter space 
where the usual waterfall condition of hybrid inflation is 
violated. 



B. Nonscale-invariant case 

In the nonscale-invariant case Jl and /nl depend on 
k. Because the tachyon spectrum is blue in this case the 
strongest constraint comes from evaluating Jnl at 
the largest values of k which are measured by the CMB. 
In deriving our constraints we conservatively take this 



to be k = e°He 



-N, 



where N e is the total number of 



c- foldings of inflation ([3]). The resulting constraints in 
the plane of log 10 g and log 10 A are shown in the right- 
hand region of figure [TJ for several values of log 10 v/M p . 
We find that the most stringent constraints come from 
/l rather than Jnl, so we expect to see distortions of 
the spectrum rather than nongaussianity at the left-hand 
boundaries of the excluded regions. (The other bound- 
aries are unexcluded for the different reasons described 
in the next paragraph.) In comparing these excluded re- 



gions to figure 12 of [15j, one sees that they are smaller 
than in our previous work, in the upper right-hand cor- 
ner. This is due to correcting an error in [15| . in which we 
failed to apply the condition (fT5)) restricting the validity 
of our analysis. 



require gv/M p < 10~ 5 ; the false vacuum energy density 
Aw 4 /4 dominates during inflation, leading to the bound 
g > 460A (v/Mp) 3 ; the reheat temperature exceeds 100 
GeV, so that baryogenesis can occur at least during the 
electroweak phase transition, leading to the lower bounds 
on A. 



C. Scale-invariant case; the adiabatic 
approximation 

On the left-hand side of figure [TJ we display new con- 
straints for which the spectrum and bispectrum are scale- 
invariant. In contrast to the right-hand side, f^L pro- 
vides the dominant constraint here, so that nongaussian- 
ity is playing the important role. To obtain these results, 
we employed a different approximation for the tachyon 
mode functions, namely the adiabatic approximation, de- 
scribed in appendix F of [ill ]. Because the tachyon has 
a small mass during the entirety of inflation (subsequent 
to horizon crossing), its mass is changing slowly, and we 
can use the standard mode functions for light fields, but 
with a time-dependent mass: 



5a(x) — 



d 3 k H 



:(-k T )^^e lkx a k + h.c. (18) 



Here r\ G = M 2 V. a a/V is the slow-roll parameter for the 
tachyon, given by 



M 2 

77, (r) = 8;?-f ln|ifr| 



(19) 




FIG. 1: Excluded regions of the hybrid inflation parameter 
space, for log 10 v/M p = — 1, —3, —5, —7, —9. 

In computing //,, /jvx over a wide range of <?, A, v/M p 
we also checked that the additional assumptions were 
respected: the tachyon mass-squared m 2 , varies linearly 
with the number of e-foldings, which was shown in 15] to 



where rj = M 2 V VV /V . We have also verified that the 
solution (fTS)) can be reproduced by flS]) in the appropriate 
limit. 

Since the mode functions have a simple form in the 
adiabatic approximation, it is possible to go farther ana- 
lytically in this case. Notably, we could find an implicit 
equation for iV* after evaluating the integral 



v/Mp 



15000 N e g 



■In 



1 + 2 x i0 6 JV, fl ^y 



(20) 



This expression for ./V* is much easier to evaluate than the 
one which arises in the WKB approximation since the lat- 
ter leads to a numerical integral f(N; A, g, v) which must 
be inverted to find the TV, which satisfies /(iV*; A, g, v) — 
v 2 /4. 

Moreover, the time (N) integral in (j4]) can be evalu- 
ated explicitly using the saddle point method, since it 
is dominated by the exponential growth near N = N*. 
Namely, an integral of the form J dNe 9 is approximated 
by e 3 " I \[\g^\ where g* is the maximum value (at N = 
iV#) and g[ is the derivative evaluated at the same point. 
Defining rjf to be \rj a \ at N = iV», this results in the 
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expression 



B(k t ) = AH^dl J ■j^\P\~ 3 ~ 2nf \P- k * 

(\ P + k 2 \- 3 - 2 ^ + b + fc!|- 3 - 2 "/) 



-3-27)/ 



(21) 



for the bispectrum, 4 where <i* = H 2 K 2 7]fe 2vf N ' /(2e). In 
the limit of small c ~ r]f, this is manifestly nearly scale 

invariant, B ~ 1/fc? + , by power-counting the diver- 
gent behavior of the integral in the infrared. The diver- 
gence must be cut off in the usual way, by ignoring modes 
with p smaller than the horizon. Numerically evaluating 
the remaining p integral for foj ~ k, and using fcr* ~ e Ni 
for modes near the horizon, 5 we find 



B(k) ^ 45 fc~ 



-6(1+77/) ( n 2 H 2 r] f e 



2r?/JV e 



27T6 



(22) 



Further using the COBE normalization to write 
K 2 H 2 /e — 2 x 1CP 7 , we find that the nonlinearity pa- 
rameter is 



NL 



-2.6 x 10" 



(23) 



Moreover the COBE normalization implies rjf = 
7360N*gM p /v, Demanding that \Jnl\ < 100 gives the 
new excluded regions on the left-hand side where is 
the mode function for the fluctuations (discussed below). 
This happens at some time after the onset of the instabil- 
ity. For a wide range of parameters (including the values 
originally considered in [3]) one has N+ <C 1 so that 
the symmetry breaking completes on a time scale short 
compared to the Hubble time. This is the usual water- 
fall condition regime of hybrid inflation. In the present 
work we will consider both the possibilities that iV* <C 1 
an also N+ > 1. of figure [TJ Unlike the nonscale- invariant 
regions, these have nongaussianity being the dominant 
effect, rather than the tachyon contribution to the spec- 
trum. 

We have claimed that in the scale-invariant regime the 
dominant constraint is coming from f^L and not from fi,, 
contrarily to the nonscale-invariant regime. We now jus- 
tify this claim. Repeating the steps above for the tachyon 
spectrum, S, one obtains 



(24) 



Thus, in the scale-invariant regime, the linearity and non- 
linearity parameters are related as 



\f NL \ ~ io e \u\ 3/2 



so that > except when \fi\ is extremely small. 

This demonstrates that it is indeed possible to obtain 
significant nongaussianity in this region of the parameter 
space. We have verified that the result (|25[) can also be 
derived using the mode function solutions (J9]). 



IV. IMPLICATIONS FOR BRANE-ANTIBRANE 
INFLATION 

We now apply our results to a popular model of in- 
flation from string theory, brane-antibrane inflation, cor- 
recting and extending our preliminary results in this di- 
rection in [l5f . This can be done by mapping the low- 
energy effective action for the brane-antibrane system 
onto the hybrid inflation model. We focus on the popu- 
lar KKLMMT scenario [H [H, [39[ which reconciles brane 
inflation with modulus stabilization using warped geome- 
tries with background fluxes for type IIB string theory 
vacua [40j. In this model, the antibrane is at the bot- 
tom of a Klebanov-Strassler (KS) throat [4l[, with warp 
factor o, < 1, and the brane moves down the throat. 6 
Within the KS throat the geometry is well approximated 
by 



a{y) 2 g ilv dx p 'dx 1 ' + dy 2 + y 2 d£l\ 



where y is the distance along the throat, a(y) = e v is 
the warp factor and dfl 2 is the metric on the base space 
of the corresponding conifold sigularity of the underlying 
Calabi-Yau space. In the subsequent analysis we ignore 
the base space and treat the geometry as AdS$. 

In the following we compute only the nongaussianity 
which is due to the preheating dynamics an d ig nore the 
possible effects of the inflaton sound speed [13, [ll|, [43| ; 
hence our results may be thought of as a lower bound on 
the nongaussianity from brane inflation. 

In string theory the open string tachyon T between a 
D3-brane and antibrane, 7 separated by a distance y, is 
described by the action j44| 

Stac = - / d A x y/^g£ 



C = V(T,y) + (a i M s )- 2 g^d IJ .T*d v T (26) 
Here the small- |T| expansion of the potential is 

V(T,y) = 

(27) 



(25) 




4 In [l5f | , the conformal time when the instability starts is (perhaps 
confusingly named) r, = —1/H, due to our choice of N = for 
the beginning of the instability. 

5 The horizon-crossing condition is kri = 1, and r*/ri = e Ni /e°. 



6 See [S for other discussions on nongaussianity in string theory 
models of inflation. 

7 We restrict ourselves to inflation models driven by D3-branes 
since inflation driven by higher dimensional branes have prob- 
lems with overclosure of the universe by defect formation |4:I . 
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where M s is the string mass scale, T3 = g~ l J (2ii) 3 
is the D3-brane tension, and g s is the string coupling. 
Notice that in the warped throat scenario the instabil- 
ity does not set in until the branes are separated by the 
unwarped string length, 8 (aiM s ) _1 . An interesting dif- 
ference between the string tachyon and that of ordinary 
hybrid inflation is that (at y = 0) the tachyon potential 
V(|T|) in the string case does not have a local minimum; 
rather 



inflation potential (TTJ). This gives the correspondence 



V(T,0)~r 3 e- 



|T| 2 /4 



(28) 



The potential is minimized as T — > 00. Therefore T 
does not have a VEV. Nevertheless, the unstable brane- 
antibrane system decays into closed strings soon after the 
instability begins, and the large- T part of the potential 
is not meaningful for determining the actual evolution of 
the tachyon. In hybrid inflation, it is also true that the 
end of inflation occurs somewhat before the fluctuations 
of the tachyon become as large as the VEV. We will see 
that even in the absence of a T 4 coupling, we can still 
define the equivalent of A and v for the brane-antibrane 
system, by equating jXv 4 to the false vacuum energy, 
and Xv 2 to the tachyon mass scale. This amounts to 
replacing the condition for the end of inflation $2$ by 



d 3 k 

W) : 



N=N* 



false vacuum energy 
jtachyon mass| 2 



(29) 



Despite that fact that the tachyon potential is only mini- 
mized at \T\ — > 00 the condition (f29|) is quite reasonable. 
Detailed numerical simulations of the symmetry break- 
ing in the theory ([26]) were performed in (45[. Compar- 
ing to the analysis of [45[ one finds that N+ as defined in 
(|29p roughly corresponds to the time at which singulari- 
ties in the spatial gradients of the tachyon field form [Z(| . 
The appearance of singularities within a finite time corre- 
sponds to the formation of lower dimensional branes [47| 
and hence by N = V* the inflaton field ceases to exist as 
a physical degree of freedom. This means that, as in our 
previous analysis, for N > iV+ there no longer exists any 
nonadiabatic pressure (since only one field, the tachyon, 
is dynamical) and the large scale curvature perturbation 
becomes conserved to all orders in perturbation theory 

The effective values of the couplings can be found by 
rewriting the action in terms of the canonically normal- 
ized fields a = ajy/j^T /M x and ip — \Jr^y (see equations 
3.6, 3.10 or C.l in [23|]), and then matching to the hybrid 



There is some confusion on the literature on this point, with 
some papers having stated that the instability is determined by 
the warped string length, but this is not the case . We thank 
L. Leblond for pointing this out to us. 



2 aiM s 



/9s 



g = a/ 2ng s m 



(30) 

(31) 
(32) 



For the analysis of the preceding sections to be valid, 
the inflaton potential must be well-described by Vb + 
\m? ip kp 1 during the relevant e- foldings of inflation. The 
full potential can be written as 



T/ m <P 2 



An 2 (p 4 



(33) 



where Vq = 2afT 3 and v is a geometrical factor which is 
given v — 27 /16 for the KS throat. It is typical to pa- 
rameterize the inflaton mass in terms of the dimensionless 
quantity (3 as m 2 v = /3H$ where H$ = V /(3M 2 ), Using 
the COBE normalization, we find that 



= l V 2 a 3 



(34) 



Demanding that the mass term in (|33[) dominate over 
the Coulomb term even when the branes are separated 
by the local string length yields a lower bound on /?: 

> 324«W 

The parameter (3 is also bounded from above by the 
requirement that \n — 1| < 10" 1 which corresponds to 
g 2 s af{M p /M s ) 2 «5x lO" 6 . 

Our results apply only in the case that /3 > 0; moreover 
the case where (3 < does not make sense from the string 
theory point of view, since ip = denotes the bottom of 
the throat, and the brane must roll toward that point, 
not away from it [49l |. 

Our preliminary results about this in [IBj suffered from 
the neglect of the condition (|13p : moreover we unduly re- 
stricted the full string parameter space by assuming that 
the scale of inflation was determined by the COBE nor- 
malization; that is not the case. As in hybrid inflation, 
we still have three parameters even after normalizing the 
spectrum, which we can take to be g s , ai and the ratio 
of the warped string scale to the Planck scale, diM s /M p . 
Taking into account the additional restrictions on /3, we 
find that the scale-noninvariant exclusions (right-hand 
side of fig. [J) do not survive at all in the KKLMMT 
model; however all the scale-invariant ones do. Therefore 
this model has the potential for producing large nongaus- 
sianity, and is even constrained by producing too high 
levels of nongaussianity. 

The constraints in the string parameter space are 
shown in figure [5] The excluded regions shown here 
correspond to very small values of g s <, 10~ 10 . In the 
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a 




-20 

lo Sio § s 



FIG. 2: Excluded regions of the KKLMMT brane-antibrane 
inflation parameter space, in the plane of log 10 at versus 



log 10 g s for log 10 aiM a /M p = 
sianity. 



13,-11, 



-5 from nongaus- 



sider the complications which arise when these models 
are coupled to moduli fields 

As in the models previously considered the false vac- 
uum energy dominates during inflation and the Hubble 
scale is given by 



H 



9t 



V6M P 



(36) 



The inflaton couples to two scalar fields of which 
one linear combination a is tachyonic. Its mass-squared 
is given by 



2/2 

9 VP 







(37) 



By comparing (|36p and (|37[) to the hybrid inflation po- 
tential |T]) we can determine the hybrid inflation model 
parameters as 



A 



£_ 

2 



2£ 



(38) 
(39) 



simplest way of connecting type IIB string theory to 
low-energy phenomenology the gauge couplings of the 
Standard Model are related to g s by running down from 
the string scale, which would render such small values 
of g s incompatible with observations. However, type IIB 
strings are dual to themselves under SL(2,Z) transforma- 
tions which take g s — > l/g s - In the dual picture, the 
string coupling is very large, and the gauge dynamics at 
the string scale would be confining. It is conceivable that 
the Standard Model arises as a remnant of a strongly cou- 
pled gauge theory at the high scale, similar to technicolor 
models. In this case, the small values of g s which give 
rise to large nongaussianity could still be compatible with 
particle physics constraints. 



The coupling g retains its original meaning in P-term 
inflation. 



A. D-term Inflation 

D-term inflation corresponds to taking / = in ([35 
During a slow roll phase the inflaton field evolves as 



g 3 C 



<t-t c ) 



which implies that m 2 varies linearly with the number of 
c-foldings. Scales relevant for the CMB left the horizon 
when ip — ip n where 



V. IMPLICATIONS FOR P-TERM INFLATION 

In realistic models of hybrid inflation from supergrav- 
ity, the potential is generated by F- or D-terms. P-term 
inflation is a class of models which combines the two 
kinds of terms and can interpolate between them [2~4| . 
The potential for P-term inflation, along the inflationary 
trajectory, is 



V int = 



1 H TT In — 77 H T 

8tt 2 tpl 8 M p 4 



(35) 



where • • • denotes terms of order ip 6 /M® and higher. 

Here (p c = £ = VK+F+W ^ s defined in terms of two 
Fayet-Iliopoulos parameters £+ and £3, and in ([55]) / 
must lie in the interval < / < 1, since it is defined 
as / = l£±J 2 /£ 2 - The limits / = and 1 correspond to 
D-term [5(| and F-term [51( models, respectively. We will 
consider each of these limits separately. We do not con- 



27r 2 M P ? 



9——M 2 

2tt 2 p 



Two distinct regimes are possible depending on the value 
of the coupling g. It is often assumed that g is rel- 
atively large so that g 2 N/(2n 2 ) > £/M 2 H|| which 
gives the correct amplitude of density perturbations with 
£ = 10~ 5 Mp and requires g > 2x 10~ 3 for consistency. In 
this regime cpN ^S> f c so that slow roll at the onset of the 
instability is not guaranteed and our previous analysis of 
the tachyon mode functions does not apply. However, in 
this regime it is also difficult to satisfy the constraints 
coming from the cosmic string tension, to avoid overpro- 
duction of cosmic strings, £ ;$ 4x 10~ 7 54]. (Ref. [55| has 
pointed out that the constraints on the cosmic string ten- 
sion can be weakened by incorporating the effect which 
strings have on the observed spectral index.) 

We are therefore driven to consider D-term inflation 
in the regime of small coupling g 2 so that g 2 N/(2ir 2 ) <§; 
i/M 2 and i^jv — fc- In this case we are guaranteed that 
the universe will still be in a slow roll phase at the onset 
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of the instability and our previous analysis holds without 
modification. This corresponds to very small couplings 
j<2x 10~ 3 , however, there is no obstruction to taking 
such a small coupling since g 2 is not necessarily related to 
the gauge coupling constant in a GUT [24| . In this regime 
the COBE normalization fixes ^~7x 10~ 4 g 2 / 3 M 2 so 
that g is the only independent model parameter. The 
cosmic string constraint (<4x 10~ 7 M 2 then restricts 
the coupling g to be smaller than g < 1.3 x 10 -5 . 

Applying our previous analysis of hybrid inflation to 
D-term inflation, 9 including the additional constraints 
mentioned above, we find that there is a range of cou- 
plings, 

- 10.0 < log 10 g < -8.7 (40) 

which is ruled out because of the spectral distortion con- 
straint, in the scale-noninvariant region of figure [TJ On 
the other hand there is no constraint coming from non- 
gaussianity for this model. 

B. F-term Inflation 



For more general P-term models with < / < 1 we 
expect the excluded regions to interpolate between ([4*0]) 
and (|4ip . In deriving our constraints we have been driven 
to the small coupling regime by the requirement that the 
cosmic sting tension be within observational bounds. Our 
analysis does not give any significant constraint on the 
string theoretic D3/D7 model [59| since in this case the 
cosmic strings are not stable and there is no motivation 
to consider the small values of the coupling g in (|40|) 
and (|4"Tj) . Indeed, such small couplings are difficult to 
motivate from string theory [60| . 

VI. THE CASE OF A COMPLEX TACHYON 

In the preceeding sections we have applied the results 
of 15], which were derived under the assumption that a 
is a real field, to models in which the tachyon is actually 
complex. In doing so we have assumed that the gener- 
alization of the analysis of [l5| to the case of a complex 
tachyon does not significantly modify the exclusion plot, 
figure [TJ Here we verify this claim. 



F-term inflation [5l|, [57| corresponds to taking / = 1 
in (|35p . In this case the dynamics are somewhat more 
complicated than the D-term model. Again there are 
two possible regimes: a large coupling regime where 
ifN 3> ifc and our previous analysis does not apply 
and also a small coupling regime where ip^ = ip c and 
our previous analysis does apply. The large coupling 
regime corresponds to g £ 2 x 10~ 3 and again the cos- 
mic strin g te nsion constraints are difficult to satisfy (see, 
however, [53]). We are therefore driven to consider only 
the small coupling regime, g <, 2 x 10~ 3 . For couplings 
3 x 10~ 7 < g < 2 x 10~ 3 it can be shown that the 
quadratic term ip 4 /M 4 in the potential (|35[) can be ne- 
glected and the dynamics is identical to D-term inflation, 
which we have already considered. Thus we consider only 
the F-term model for g<3x 10~ 7 since this is the only 
region of parameter space for which the model differs sig- 
nificantly from the D-term model. 10 

For g <, 3 x 10~ 7 , so the /-term is dominating the po- 
tential, the slow roll parameter e = £ 3 /(8Mp), and the 
COBE normalization fixes £ = 6.7 x 10 6 g 2 M 2 and the 
cosmic string tension is within observational bounds for 
g <, 2 x 10~ 7 . Again applying the general hybrid inflation 
constraints, we find the excluded region 

-13.0 < log 10 .g < -9.5 (41) 

which, as in the case of D-term inflation, comes from the 
A; 3 spectral distortion effect rather than nongaussianity. 



See f56tl for further discussion of preheating in D-term inflation. 
We have neglected the intermediate regime 0.06 ,§ g ^ 0.15 
which will not yield significant nongaussianity or spectral dis- 
tortion. 



A. Cosmological Perturbation Theory for an O(M) 
Multiplet 

Before restricting to the case of a complex tachyon we 
consider the somewhat more general case of an O(M) 
symmetric multiplet of tachyon fields a a with A = 
1, . . . , M. The matter sector is expanded in perturba- 
tion theory as 

<p{T,2) = ^ (0) (r)+5 (1 V(^^ + ^ (2 V(^^(42) 
a A {r,x) = 6^a A {T,x) + ^a A (T,x). (43) 

As in [HI the time-dependent vacuum expectation value 
(VEV) of the tachyon fields are set to zero {a a) = = 
which is a consequence of the O(M) symmetry of the 
theory. Notice, however, that the tachyon field does de- 
velop an effective VEV for the radial component 

(H) = (V^ 1 ) 

We also assume that 

dV _ d 2 V _ 
da a d/JAdtp 

but V is, for the time being, otherwise arbitrary. Here 
and elsewhere the potential and its derivatives are under- 
stood to be evaluated on background values of the fields 
so that V = V(tp(°\ cr (0) ), for example. 

We consider only the 5 (2 ^G% = n 2 5^T^, d l 5 { ^G\ ) = 
K 2 diS^T^ and 8)5^Gl = K^SjS^T? equations since 
the second order vector and tensor fluctuations decouple 
from this system. In the case that er^ = the second 
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order tachyon fluctations 5^ a a decouple from the infla- 
ton and gravitational fluctuations up to second order and 
hence we do not need to solve for 5^ a a- Note also that 
the Klein-Gordon equation for the inflaton fluctations is 
not necessary to close the system. In this section we 
sometimes insert the slow roll parameters e and r\ explic- 
itly though we do not yet assume that they are small. We 
also introduce the shorthand notation m 2 = d 2 V/dip 2 . 
The second order (0, 0) equation is 

37^' (2) + (3 - €)H 2 ^ _ d k d k ^ 

Ti (44) 



K 

y 



dip 



where Ti is constructed entirely from first order quanti- 
ties. Dividing Ti into inflaton and tachyon contributions 
we have 

Ti = Tf + TI 

where 

Tf = 4(3-e)H 2 (0«) 2 + 2 K Vo0 (1) <5 (1 V' 

V 



K 2 / \ 2 



\a 2 m% - K - (v^V) 

+ 8^d k d k ^ +3^Y + 3(V^) 2 (45) 



and 



l>2 V -S^a A S^a B 



da Ada B 

The divergence of the second order (0, i) equation is 



(46) 



d k d k 



#W + n^A = ^d k d k 5™ip + T 2 (47) 



where T 2 = T 2 + T 2 is constructed entirely from first 
order quantities. The inflaton part is 



and the tachyon part is 

The trace of the second order equation is 



(48) 



(49) 



3^" (2) + d k d k 

m<t>' {2 

3^ " 



<A (2) - 1/> 



(2) 



6Hip 



/(2) 



3(3 - e)H 2 (f> (2) 
9 



(50) 



where T3 = T3 + T3 is constructed entirely from first 
order quantities. The inflaton part is 



12(3 - e)H 2 (V«) 2 - 6kVo0 (1) * (1 V 



3k fdi) A 2 3 « 2 

+ — ( 5i y) - t 

+ 3 (V«) 2 + 8<^9 fe 9^« + 24W^V (1) 
7(W«" 2 
and the tachyon part is 



(51) 



^a' A 5^a' A - \d l 8^a A d l 8^a A 
, 8 2 V 



2 d<7Ad<JB 



5 w a A 8 {1) (JB 



(52) 



The derivation of the master equation which was pre- 
sented in appendix B of [15] follows here unmodified ex- 
cept for the new definitions of Tf , T 2 and Tf. The 
master equation is 



4>"^ + 2(7? - e)Hcb' {2) + [2(r] - 2e)H 2 - d k d k ] 
where the source is 



= J 
(53) 



J = Ti-T 3 + 4A- 1 T 2 + 2(l-£ + r7)HA- 1 T 2 
+ A-y - (1 + 2e - 2r))HA- 1 -/'. (54) 



and the quantity 7 is defined as 



7 = T 3 - SA^Ta - 6HA^T 2 



(55) 



We can split the source into tachyon and inflaton contri- 
butions J = J v + J a in the obvious manner, by taking 
the tachyon and inflaton parts of Ti, T 2 , T3, 7. 

In appendix B we prove the identity (see eqn. IB-2[) 

7. = -y (d^a A d^a A ) 

- 3k 2 A- 1 ^ (d k d k 5^a A d l 8^a A ) 

which is analogous to the result for a real tachyon field, 
derived in (la ). 

We now proceed to derive the tachyon curvature per- 
turbation. The derivation of presented in [l5| follows 
unmodified except, of course, for the change in the defi- 
nitions of T^, T 2 , T3 and 7^. From this point onwards 
we assume that e,\rj\ <C 1. The leading contribution to 
the tachyon curvature perturbation is 



A2) ~ 



dr' 



T? 



1 Tf 



H{t') 3H(t') 



2H(t'Y 



3 H{t) 



3" 1 3 
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Now, using equations (|46|) and (|52f we can write this in 
terms of the tachyon fluctuation S^'a as 



Ca> = - I dr' 

e J-l/aiH 



5^a' A S^a ,A 



H{t) 3 
d 2 V 



da A da B 



H{t>) 
5 {1) a A 8^<j B 



(56) 



The corrections to ([55)1 are either total gradients or are 
subleading in the slow roll expansion. In deriving (|56|) we 
have restricted ourselves to the preheating phase during 
which the fluctuations 5^a A grow exponentially. 

Using ([56]) the second order tachyon curavture pertur- 
bation can be computed once the fluctuations 5^<ja & re 
determined. The first order tachyon fluctuations are de- 
scribed by the perturbed Klein-Gordon equation 



OA ■ 



d k d k S {1 



, d 2 V 
d(JAdo B 



5 W 



<?B=0 

(57) 



B. Complex Tachyon Mode Functions 

At this point we restrict our attention to the case with 
M = 2 and the potential 



V 



9 2 2 A 



(58) 



For a 



A — the mass matrix is diagonal 
d 2 V 



da Ad a b 



= (-Aw 2 + g 2 tpl) 5 A b 
= m 2 5 A B 



so that the tachyon fluctuations with A = 1 and A = 2 
evolve independently (see eqn. [57)) . 

As previously the quantum mechanical solutions 
5^'ga are written in terms of annihilation and creation 



operators a A , a, ^ in the usual way 



5 {1) <ja(x) = 



d 3 k 
(2tt) 3 /' 



a£Ut) e ikx +h.c. 



(59) 



Both components A — 1 and A = 2 have the same time 
dependence owing to the fact that the mass matrix is di- 
agonal. The £fc in (f59"|) are thus identical to the solutions 
of ([H]), which we have already studied. 



C. The End of Symmetry Breaking 

For a multi-component tachyon the condition defining 
2V* must be modified as (S ( ^a A 8 W (T A )(N = N*) = v 2 /A 



which, for the case M — 2, changes ([2]) to 



d 6 k 
(2^)3 



N=N i 



D. Tachyon Curvature Perturbation 



For the potential (f58"| the tachyon curvature perturba- 



tion (i 2) 



decomposes into a sum of term 

A 2 ) = 



E 

A=l,2 



A 2 ) 



where C^ A is the contribution to 



Consider, as an example, 
curvature perturbation 

-(2) A 2 ) 



coming from a a- 
the spectrum of the tachyon 



//■(*) m z j \ _ i A 2 ) A 2 ) \ I /A z > /-( z > \ 



42) A 2 ) 



,A2) A 2 ) \ , (sw c w \ 

Because the annihilation/creation operators a\ and a 2 
are independent the cross-terms on the last line do not 
contribute to the connected part of the correlation func- 
tion. This means that 

/A 2 ) A 2 ) \ _ o /A 2 ) A 2 ) \ 



(2) A 2 ) 



The quantity (( {2) kl {[%) = (C^L^O wiU be identical 

(2) (2) 

to the (Q k Q k 2 ) which we have already computed. We 
see, then, that the effect of having a complex tachyon 
field (as opposed to a real field) is to multiply /l and 
Jnl by a factor of 2 and also to slightly reduce . The 
net change in /j,, Jnl is order unity and the new exclu- 
sion plots is difficult to visually distinguish from figure 
[TJ This justifies our previous claims that our constraints 
do not change significantly when one considers a complex 
tachyon field. 



42) ,(2) 



VII. CONCLUSIONS 

In this paper we have studied the evolution of the sec- 
ond order curvature perturbation during tachyonic pre- 
heating at the end of hybrid inflation. We have found 
that, depending on the values of certain model parame- 
ters, two interesting effects are possible: 

• Preheating generates a scale-invariant contribution 
to the curvature perturbation. In this case signif- 
icant nongaussianity can be generated during pre- 
heating and the model is even constrained by pro- 
ducing too high a level of nongaussianity. 

• Preheating generates a nonscale-invariant contribu- 
tion to the curvature perturbation with spectral in- 
dex n = 4. In this case the strongest constraint 
comes from the distortion of the power spectrum 
and no significant nongaussianity can be produced. 
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In both cases one typically requires fairly small val- 
ues of the dimensionless couplings g, A in order to ob- 
tain a strong effect. Note that a small coupling g 
does not require fine tuning in the technical sense, 
since g 2 is only multiplicatively renormalized: f3(g 2 ) ~ 
0(g 2 \, g 4 )/(167r 2 ). That is, if g is small at tree level 
then loop corrections do not change its effective value 
significantly. 

We have applied our constraints on hybrid inflation to 
several popular models: brane inflation, D-term inflation 
and F-term inflation. In the case of brane inflation we 
have found that significant nongaussianity from preheat- 
ing is possible for sufficiently small values of the warp fac- 
tor. For both D- and F-term inflation we have shown that 
no nongaussianity is produced during preheating, how- 
ever, we still put interesting constraints on the model due 
to the distortion of the spectrum by nonscale-invariant 
fluctuations. 

We have also generalized the results of [15| to the 
case of a complex tachyon field, confirming our previ- 
ous claims that this modification does not significantly 
alter our exclusion plots. 

We should note that the model of hybrid inflation con- 
sidered here always gives a small blue tilt to the spectral 
index, n > 1, which is disfavoured by recent data ||. 
One avenue for future study 61[ is to generalize our re- 
sults to the case of inverted hybrid inflation [g^] which 
always gives n < 1. 
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and may be written as 

{\W-x for the branch with iV fc < -1; 

±Wo (~^) for the branch with -1 < N k < 0; 
|W (+^r) for the branch with N k > 0. 

(A-2) 

In the region k > ^Jc/(2e) the solution is single valued 
and can be written as 



(A-3) 



One may derive some asymptotic expressions for N k in 
various regions of interest. When |JVjt| > 1 we have 



(A-4) 



which describes N k at k ^> -J cj (2e) and also the lower 

branch of N k at k < ^/c/(2e). For k < >/c/(2e) there 
are two more branches of the solution with approximate 
behaviour 



k 2 

N k ~±- 



(A-5) 



In our analysis we have used the approximation that 
N k is a single- valued function, described by 




+ -e(k- V^/W)) Wo U-j 

where 0(x) is the Heaviside step function. We have 
verified both numerically [l5j and analytically that the 
single-valued approximation does not significantly alter 
our results. 



APPENDIX B: An Identity Concerning 7 CT 



APPENDIX A: The Matching Time N k 

The matching time N k which determines the bound- 
ary between large- and small-scale behaviour of the mode 
functions ([9]) is determined by the transcendental equa- 
tion 



\N k \e 



2N k 



c 



(A-l) 



The solutions may be written exactly in terms of the 
branches of the Lambert-W functions. In the region k < 
\J c/(2e) the solution is triple- valued (see figure 1 of [l5j ) 



In this appendix we derive an identity concering the 
tachyon source term 7 CT 



= Tg — 3A" 1 9 T T^ - 6WA- X T£. 
Using equations (|49|) and (|52)l we can write this 



^a 2 -^-d k d k (S^a A S^a B ) 
6Hd i {S^a' A d i 5^a A 
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and, after some algebra, we have 



l -d k d k {d % 8^a A d l 5^a A ) 



which follows from the 0(M) symmetry of the theory. 
The last two lines of (|B-lj) can be simplified using the 
^equation of motion for the tachyon fluctuation (|57|) which 



gives 



35, \5^a' A + 2H5 {1) a' A 



, d 2 V 
da Ada b 
8 2 V 



OB 



6 {1) a B d k d k S^a J 



da A da B 

In deriving this equation we have used that fact that 
d 2 V d 2 V 



da Adas da b da a 



7. = -y (d^a A d^a A ) 

- 3K 2 A- 1 d i (d k d k 6^a A d i S^a A ) (B-2) 
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